A Very Narrow Shadow Extra Z-boson at Colliders 
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We consider the phenomenological consequences of a hidden Higgs sector extending the Standard 
Model (SM), in which the "shadow Higgs" are uncharged under the SM gauge groups. We consider 
a simple (7(1) model with one Higgs singlet. One mechanism which sheds light on the shadow sector 
is the mixing between the neutral gauge boson of the SM and the additional (7(1) gauge group. The 
mixing happens through the usual mass-mixing and also kinetic-mixing, and is the only way the 
"shadow Z" couples to the SM. We study in detail modifications to the electroweak precision tests 
(EWPTs) that the presence of such a shadow sector would bring, which in turn provide constraints 
on the kinetic-mixing parameter, s e , left free in our model. The shadow Z production rate at the 
LHC and ILC depends on s e . We find that observable event rate at both facilities is possible for a 
reasonable range of s e allowed by EWPTs. 



I. INTRODUCTION 

In the pursuit of physics beyond the Standard Model 
(SM) it is very common to encounter one or more Abelian 
gauge symmetry than the SM [7(1) hypercharge. Two fa- 
miliar examples are the grand Unified theories (GUTs) 
based on SO(10) that breaks to Gsm X £7(1), where Gsm 
is the SM gauge group and £'(6) which ultimately breaks 
to Gsm x U(l) x [7(1). Because of their GUTs parent- 
age the extra Z bosons from the breaking of the U(l) 
symmetries have tree level couplings to the SM particles; 
in particular the fcrmions. This makes them highly vis- 
ible and their phenomenology has been well studied 0. 
More recently extra dimensional models with extra [7(1) 
gauge symmetries in the brane world scenario are increas- 
ingly popular. A feature of this newer construction is 
that the extra [7(1) factors can be hidden from the visi- 
ble sector. Hidden sectors are motivated also by studies 
in supersymmetry breaking mechanisms. Independent 
of the theoretical motivation, extra Z bosons from hid- 
den sector typically do not have direct couplings to the 
SM particles. Their phenomenology can be very different 
from visible extra Z's. They are also harder to produce. 
As the start up of the LHC draws near, the search for 
extra Z bosons is a high priority item due to their rela- 
tively clean signatures from Drell-Yan processes. Clearly 
it is important to include hard to find extra Z bosons 
in this search. Although these bosons have no direct 
couplings to SM particles, they can manifest themselves 
through mixings with the SM Z boson, and so are not 
completely invisible. Since the mixing is crucial for phe- 
nomenology we construct the simplest model of this kind 
to capture the physics of such an extra Z boson. It has 
the gauge symmetry Gsm X U(l) s where the subscript 
denotes "shadow" : the name will become clear later. The 
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SM fermions are singlet under U(l) s . This [7(1) s is bro- 
ken by a shadow Higgs sector which is just the Abelian 
Higgs model with a complex scalar <f> s . The (j) s field is 
a sinlgct under Gsm but interacts with the SM Higgs 
bosons via renormalizable interactions. The complete 
Lagrangian is given by 



C = C 



SM 



-X^X„ 



2 M 



(1) 



where B^ v is the field strength tensor of the SM hyper- 
charge U(1)y, $ is the SM Higgs field, and g s is the 
gauge coupling of the shadow U(l) s - For simplicity we 
have normalized the shadow charge of <fi s to unity. The 
kinematic mixing of the two [7(l)'s is parameterized by 
e, which a priori need not be a small number. For a 
visible extra Z this mixing term is expected to be only 
induced at the loop level [2j, and thus |e| <C 1 is gener- 
ally assumed in its phenomenological studies How- 
ever, this need not be the case here. Indeed, a calculation 
in string theory of the mixing-generating vacuum polar- 
ization diagram shows that in general, one can expect 
kinetic mixing effects on the order oflO -4 ~ 1CP 2 at 
the weak scale (barring accidental cancellations in the 
tree level spectrum) jjj. Given the theoretical signifi- 
cance outlined above, we shall leave e as a free parameter 
to be constrained by experiments in particular the elec- 
troweak precision tests (EWPTs). Now it is well known 
that the kinetic terms including the mixing can be recast 
into canonical form through a GL(2) transformation. Ex- 
plicitly, this is given by 



c e 

-8 C 1 



where 



X' 
B' 



1 



(2) 



(3) 



After spontaneous symmetry breaking (SSB) X' and B' 
will mix resulting in a shift in the SM Z mass. The 
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physical bosons are now linear combinations of the two. 
The photon will remain massless, and the W bosons will 
be unchanged from the SM. This is expected since the 
shadow sector only interacts with U(l)y through the 
shadow Higgs interactions. The details of this symmetry 
breaking is given in Sec.[H] Feynman rules for the model 
are also given there. In this paper we focus on the phe- 
nomenology of the physical shadow sector neutral boson, 
Z s . Since we are interested in collider physics we shall 
assume that vacuum expectation value (VEV) of (f> 8 is 
of the order of the weak scale or higher. In Sec. II I II we 
study the impact Z s has on electroweak precision mea- 
surements. From these, as well as anomalous magnetic 
moment of the muon and recent results from M0ller scat- 
tering, we derive constraints on the parameters of our 
model, in particular on s e . We employ a conservative 
strategy and demand that the fits to the data are not 
much worse than that of the SM. With these limits in 
hand we explore the prospect of observing the Z s at the 
LHC and the ILC in Sec. lIVI Finally we give our conclu- 
sions in Sec. [V] Recent work with an extra Z similar to 
ours is given m and the older literature can be found 
in 



SYMMETRY BREAKING AND THE 
SHADOW WORLD 



rotation: 



II. 



The most general renormalizable Gsm X U(l) s invari- 
ant scalar potential is: 

+/. t 2 $ t $ + A($ t $) 2 . (4) 

This Higgs potential is also used in phantom Higgs mod- 
els [U- After SSB the scalars acquire nonzero VEV, 



V2 



with 



V 2 



AA S 



(5) 



(6) 



To ensure that the potential is bounded from below and 
the above values correspond to a global minimum we re- 
quire A, A s > and k > 0. 

The SU(2)lX U(1)y x U(l) a symmetry is broken down 
to U(1)qed- This pattern of breaking is peculiar in 
that the mass of the W boson remains as in the SM, 
i.e. M w = (g 2 v )/2. In the neutral sector we have a 
massless photon and two massive neutral bosons which 
are not yet in the mass eigenbasis. The usual SM defini- 
tion: tan 6*vi/ = 3yIqii electric charge e = 52 sin 6w, an d 
Q^ L r = R + Y[ R remain intact. 

For the neutral gauge bosons the transformation be- 
tween the weak and mass basis is given by the following 
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S ,j 


c v J 




\ z s 



where sw (cw) denotes sinOw (cosByy) and similarly for 
the rotation angle 77. The first rotation is the standard 
one that gives rise to the SM Z and the second one di- 
agonalizes the mixing of the two Z bosons. The mixing 
angle is given by 



tan 2rj = 



2sws e 



c 2 w (M 3 /M w ) 2 



1 " 



(8) 



where M 3 = (g s v s )/2. For small s e and cwM 3 > M\y, 
rj < e. The masses for the two massive neutral gauge 
bosons are readily found to be 



M 



z,z s 



AT 2 
2c 2 



f c w M 3 



1 



'We 



yf [(cwMs/Mw) 



1 



(9) 



For the case where M 3 > Mw the Z-Z s mixing is pro- 
portional to s c , which is related to the amplitude of the 
kinetic mixing term B^X^. 

The most stringent constraints on any extra Z model 
come from EWPTs, and so we consider next the gauge 
fermion couplings. These can be readily read off from 
the Lagrangian. For the photon (A^), the SM result is 
retained as it should: 



A»Jf : i^eQ f . 



(10) 



For Z, Z s , the coupling are slightly different from the 
SM, but still flavor universal: 



Z"ff ■ H^— \{c v gf - s v s w s e Y f L ) L 

Cw ^ 

+ (c n gf - s n s w SeY f R ) R 

Z?ff ■ \(s n 9f ~ c v s w s e Y f L ) L 

cw L v J 1 



s n gf - c n s w s e Y f R ) R 



(11) 



(12) 



where g f L R — T 3 (fL t n) — s w Q f is the coupling of the 
SM Z to fermions. We see that the neutral current cou- 
plings are not only rotated as indicated by the c v fac- 
tor, but also contain an extra piece proportional to the 
fermion hypercharge due to U(l)-U(l) s mixing. Hence 
we need to reexamine the electroweak precision data us- 
ing the full couplings as well as taking into account the 
effects due to virtual Z s exchanges. For the pure gauge 
sector, it's straightforward to work out the Feynman 
rules. For example, the Feynman rules for 4 gauge boson 
(V^V 2 W+W-) vertex read: 



- igl C 4 U 2 [2gnug\p - g^xgup - g^ P gv\] 

where the factor C4 are listed below: 



(13) 
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w- V? 



W+(k 2 ) 



FIG. 1: The triple and quartic gauge vertices and the mo- 
mentum labelling. 
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Similarly, the Feynman rules for triple gauge coupling for 
V^kijW^ (k 2 )W^ (k^), where all momentum fc's go into 
the vertex, read: 

~ig 2 C^ [(fci - k 2 )\g l _ iv + {k 2 ~ k s ) ll g uX + {h - ki) u g Xll ] , 

(14) 

and for different Vs the CVs are 
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Z 


Z s 


sw 


CW c r} 


— CW s r) 



Now a few remarks about the scalar sector. The SM 
Higgs doublet has 4 degrees of freedom (DOF) and the 
shadow scalar has 2 DOF. After SSB in both sectors, one 
DOF of each scalar becomes massive physical scalar. So 
we are left with 3 + 1 massless DOFs which will be eaten 
by two W, one Z, and one Z s , and the DOF budget is 
balanced. Therefore, the shadow world gives us one extra 
neutral scalar and no charged scalars. How heavy they 
are is an interesting question. We will assume here the 
lighter one is SM-like and has a mass greater than 114 
GeV and the heavier one is more than 200 GeV. This 
amounts to assuming v s > vq and no fine tuning of the 
Higgs parameters. In the basis of {h% M , h®}, the mass 
matrix for these two neutral scalars is 



Xv 2 KV V S 



I KV V S \ p V 2 I 

It can be diagonalized by a rotation 



(15) 



fermions, the scalar-fermion couplings are given by 

g 2 m f 



hiff ■ — «cosa 
h 2 ff '■ — «sina- 



V2M- 



w 



(19) 
(20) 



\/2M w ' 

and for the gauge-scalar couplings one has 

,0™™ ig^Mw , , 2 „„ , 01 .. 

h^Z^Z : g °a\ c V + SwSeSr,) 3 , ( 21 ) 

c w 

h 1 Z^Z s : 2 c a (s v - swSeCri) g** , (22) 



c w 

X(c, ; + SwScS^iSrj - SwSeC v ) , (23) 

h\W + ^W-' v : ig 2 M w c a g^ . (24) 

For h 2 , the "shadow Higgs", c a is replaced by s a in the 
above expression. 

III. PHENOMENOLOGY 

We now perform a systematic phenomenological study 
using the Feynman rules derived previously. We present 
the analytical result involving parameters of the shadow 
sector; numerical values are summarized in Table [I] at 
the end of the section. We begin with the anomalous 
magnetic moment of the muon. 

A. Muon g - 2 

The one-loop Z boson contribution to the anoma- 
lous magnetic moment of a charged lepton is now dif- 
ferent from the SM due to modification of its coupling to 
fermions. Furthermore, the same one-loop diagram with 
Z s running in the loop contributes as well. Plugging in 
the gauge-fermion interactions obtained in the previous 
section, the SM one-loop Z boson contribution now shifts 
by an amount: 



h2 



cos a sin a 
— sin a cos a 




and the mixing angle satisfies 



tan(2a) = 



2kvqv s 



X,v 2 - Xi 



with mass square 

l r. 



'1,2 



(16) 



(17) 



X s v 2 + Xv 2 ± J (X p v 2 - Xv 2 ) 2 + 4k 2 v 2 v 2 



(18) 

The Feynman rules for the scalar sector can be read- 
ily worked out. For instance, in the mass basis of the 



5a„ = slAai M 



2 2 
.92 m £ 



1 - 



Mz 
M z 



(25) 



87r 2 M 2 c 2 w 
87r 2 M 2 c 2 w 



1.94 x 10" 



s 2 
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2 2 2 
s W S e S rj 



2 1 \ 1 2 2 2 



M 2 



[s 2 ~ 0.2354a„c,* e + 0.1850s 2 s 2 ] , 



(26) 



where we have set s 2 ^ = 0.2311. The W boson and 



other SM contributions remain the same. The two Higgs 
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bosons and higher loop diagrams contribute negligibly. 
We shall see later that Sa^ above does not give better 
constraints than EWPT. 



C. M0ller Scattering at SLAC 

The SLAC E158 M0ller scattering experiment [TJ 
measures the parity violating asymmetry, 



B. NuTeV 

The NuTeV experiment measures the ratio of neutral 
current to charged current cross-sections in deep- inelastic 
i^-nucleon scattering 0. As was suggested by Paschos- 
Wolfenstein ^(J to reduce theoretical and systematic un- 
certainties, the precision observable to measure at NuTeV 
is 



Rpw = 



{9 N L ? ~ 



-» fiX) - aiy^N 
(n N \ 2 - - - s 2 

\9R ) — r, S W 



(27) 



where (g£ R ) 2 



(9l 



L,R> 



(9t 



L.RJ 



Since the shadow 
world affects only the neutral current processes, presence 
of a new neutral gauge boson will only affect the numer- 
ator of Rpw- For the v^q elastic scattering process, the 
squared amplitude receives corrections from the modifi- 
cations in the couplings of the SM Z, and from contribu- 
tions due to the exchange of the virtual Z s . Incorporating 
these effects due to the Z s , a straightforward calculation 
shows that the numerator of Rpw is proportional to 



q—u,d 



M 2 

1^-9l[<9l-x q R 9 R ] 



' M l 



-O 



(Ml 
\m* 



,(28) 



.4 



PV 



&l — or 

<JL + Or 



(30) 



at momentum transfer Q 2 = 0.026 GeV 2 . The subscripts 
L and R denote the incident electron polarization. At 
tree level, the asymmetry is, up to C(g|): 



.4 



G F s 



PV 



- y) 



V2vra 1 + y 4 + (1 - y) 



-(<?? 



9r), (31) 



where y = Q 2 / s ~ 0.6. 

The denominator in the above expression represent the 
leading M0ller cross-section due to photon exchange, and 
the numerator is the parity violation due to photon-Z in- 
terference. It is easy to extend it to include the Z s contri- 
bution. We only need to keep the photon- Z s interference 
term: 



.4 



G F s 



2/(1 - y) 



PV 



V2na 1 + y 4 + (1 - y)- 



-(*' 



c* 2 ) (Mil 



(32) 

Assuming that Mz s ^S> Mz, the Z s effect can be ig- 
nored. From the modified Zee coupling we have 



SA 



PV 



A 



PV 



This translates into 



3 s ^ s w s 1 
1 1 _ s 2 

* 4 b W 



s 2 



(33) 



nv 



5 sin 



9.171s 2 s 2 + 18.60c v s v s e ] . (34) 



As 2 w SA PV 



4 A PV 
0.019 [s 2 + 9.17s 2 s 2 + 18.60c^s,,s e ](35) 



where x{ R = —(s^g^ R + c v Y$ sws t ) is the Z s ff cou- 
pling given in Ea. ((l"2*jl . and q 2 <C M§,M§ is assumed, 
with q the momentum transfer. The first term in Eq. I|28|) 
is the SM result, while the second term comes from the 
SM-shadow interference; we have ignored the term sup- 
pressed by (Mz/Mz s ) 4 - Note that for the isoscalar tar- 
gets considered at NuTeV, the sum is over u and d quark 
distributions. Assuming that Mz s 3> Mz, only the SM 
contribution need be kept while taking into account the 
modifications to the Z-fermion coupling. In terms of the 
mixing parameters of our model, the effective nucleon 
coupling to the SM Z is given by: 

(9l f * (9l)sm [1 - H + 1.0028 We ] , 

(9r) 2 - (9r) 2 sm [1 - 24 + 5.1218 We ] . (29) 

Note that the coupling to neutrinos is absorbed into the 
effective couplings here. 



D. Atomic Parity Violation 

In the atomic system, the exchange of SM Z boson will 
generate the parity violating Ml transition. This opti- 
cal line can be accurately measured and used to compare 
with the theoretical prediction [12J . Since the momentum 
transferred by the Z boson is much smaller than nuclear 
mass, it can sense the weak charge of all the quarks co- 
herently. The relevant quantity is: 

Qw = -4(si - g%)[(2Z + N){gl + g R ) 

+ (2N + Z)(g d L +g d R )}- (36) 

The expression for the contribution from an extra neutral 
gauge boson, X, is same as above with gL,R. changed to 
the corresponding couplings for X and multiplied by an 
extra mass factor m 2 z /vn 2 x . 

If Z s is much heavier than the SM Z, its tree-level ef- 
fect goes like s 2 (Mz /Mz s ) 2 , which can be again ignored. 
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Therefore, the leading change to Qw comes from the 
modification to the couplings of SM Z boson to fermions. 
At tree level, we have for Cs^ 3 



SQ 



~ -si + 0.7605sls 2 + 2.0627c I? s„s £ , 
Qw 



and for Tl 2 f 



SQ 



ir 



AY 



(37) 



-si + 0.7195sls 2 e + l.9774c v s v s e . (38) 



E. Asymmetries at LEP 



Consider first the general expression of differential 
cross section for the process e~ + e + — > / _ + / + mediated 
by more than one neutral gauge boson. If we ignored all 
the light fermion masses, the differential cross section for 
/~ deflected from the incident e~ direction by angle 9 
are given by: 



- A "" < " KxY + cos61 ) 2 (?xl9yl9xl9yl + 9xr9yr9xr9yr 



d(cos9) 8s 



X.Y 



-(1 - cos 0) 2 (9xl9yl9 X r9yr + 9xr9yr9xl9yl) } . ( 39 ) 



where indices X, Y run over all neutral gauge bosons and the K's are kinematic factors. For instance, A" 77 = Q 2 , and 
forX,y^ 7 , 

f-2Q f \ a{s-M%) 



J W°W 



J [{s-M x f+M x T x Y 



XX ic^ w ) (s - M x f + M%T* ' (41) 
K = ( 2 \ s 2 [(s-M 2 ){s-M 2 x )+M Y T Y M x Tx] 
XY Kc^w) [(« - M Y f + M Y T Y ][(s - M\Y + M X T X ] ' [ 1 

where M and T are the mass and width of the neutral gauge boson respectively. For photon coupling, it has been 
normalized to be — 9r = 1 ■ For other neutral gauge boson coupling, the coupling is normalized by the SM strength 
{92 /cw)- The forward-backward asymmetry is then given by 

, 3 Ex,y k xy(9xl9yl ~ 9xr9yr) 9xl9yl 9xr9yr ) 

a fb = J ■ )— f -4 • ( 43 ) 

Ex,y k xy(9xl9yl + 9xr9yr) \9xl9yl + 9xr9yr) 

In SM, the fermion's left-right asymmetry can be derived from the left-right forward-backward asymmetry: 

A f =± i*-i B -i* +a r , (44) 

6 °lf + a LB + a RF + °rb 

where L (R) stands for the left(right)-handed incident electron, and F (B) stands for the forward (backward) direction, 
cos 9 > (< 0), of the final state fermion. When more than one massive neutral gauge bosons are present, the effective 
left-right asymmetry becomes: 



a J2x,y k xy(9xl9yl + 9xr9yr) \9xl9yl ~ 9 X r9yr) ^ 
J2x,y k xy(9xl9yl + 9xr9yr) (?xl9yl + 9xr9 Y r) 



At the Z pole, K Z z = 42298.1. The other K factors are 
very small and all of them can be ignored except Kxx, 
which has a very narrow and high spike when Mz s ~ Mz ■ 
However, it drops very quickly when Mz 3 fall outside 
the Z width. Therefore, for a heavy Z s we only need to 



consider the modification of the SM Zff coupling and 
its effect on the precision measurement. 

In Table [fl we summarize t he c urrent LEP, NuTeV, 
and SLAC M0ller status (from 0]) and our prediction. 
The second column A exp is the experimental fractional 
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deviation from the SM prediction and the combined the- 
oretical and experimental uncertainty, (5 exp ), is shown 
in the parenthesis. The third column gives the fractional 
deviation from the SM our shadow model predicts. 



Quantity 


A _ /KX P \ 

p \ bJVl ) 


1 As= (Model) 1 


r z 


a nnns/'i ^ 


— S„ -+- U.O / OUSt? C^Sg 


Chad 


1 n nm 7i'q'\ 
+u.uui / i^yj 


1 a no^QQ *> ,n 
+U.UzDyos T) c, 7 s e 


1 yilCLCL ) 




,-, - 1 n /ioo7o ^ 
— -\- U.4oZ / Cf}S € 




n nn^f^n 1 * 


c 2 1 n QC1 « 


1 {11) 


n nnn/is/'i n7'\ 


J_ n 7QQQ0 
— Sjj — h U. i OiJoSTjCrjSt 




_i_n nn9fi/'9Pi s \ 


a oncic « 
— U.oUDOS^C^Sg 


i ' ft 


-i_n nni fi/'i 




D 

-Or 






R, 
lit 


-i_n nn^/i f^i 


j_n n^7fie ^ 


rt c 


n nni o. (a 7 a \ 


A 1 Qflfic r- q 


4 (0,e) 


n 1 ns/'i ka \ 


— 00. D ( St? C^Sf 


4 (0,M) 


j_n n^Qfs^ 


— 00. D ( SfjCrfSe 


4 (0,r) 
FB 


in 1 cc/i a*^ 


— 00. D / Srj C^S e 


4 (0,6) 
^FB 


n n^/i f 1 7^ 




A (0,c) 
FB 




91 9/1 c r> q 


A (0,s) 
FB 


n at Pi/'l 1 1 \ 


1 q f;Qc ^ 


1 

A e 


1 A AOO ( 1 £i\ 

+U.Uz5(lbJ 


— Ly.odQSr/CrjSe 




I A A/1 Af' /I Q"\ 

+U.U4y[4z J 


1 A QQCf, „ „ 

— L\).OODS v Cr]S £ 






— ± y .OOui)j (--77 A e 
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a nor /1 no 1 
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J T 


n n7fUi no^ 

— U.U / O^IUZ ^ 


1 Q ^^^e r> c 
— iy .OOOSrj C-qSe 




— U.UZZ^OU J 


1 Q ^^^e r> q 
— iy .OOOSj) CrjS e 


A, 

* *-b 


U.U1UI ^1 J 


— 9^1 1 r Q, 

\J . £l\J _L (5 77 O £ 


A c 


+0.003(39) 


— 1.909s^c^s E 


A s 


-0.043(97) 


— 0.251s,, CjjSe 


NuTeV 






(9l) 


-0.013(5) 


-2s 2 + 1.0028s^c^s E 


(9r) 


+0.023(36) 


—2s 2 + 5.1218s^c^s e 


SLAC M0ller 






sin 2 Oeff 


+0.007(22) 


0.019s 2 + 0.353s^c^s e 


Qw(Cs 133 ) 


-0.0068(66) 


-s 2 + 2.0627c„s„s e 


Qw(Tl 205 ) 


-0.0018(317) 


-s 2 + 1.9774c^s^s e 



TABLE I: The comparison of experimental values and the- 
oretical prediction, up to C(s 2 ), for various EWPT observ- 
ables. 

To give a measure of how well our model fits the data 
from EWPTs compare to that using purely the SM, we 
define a number, that measures the deviation between 
theory and experiment, 



A 



exp 



J exp 



For SM only, the deviation is 



xf f = X2(0,Af 3 Hx 2 (s £ ,oo) 



A! 



34.908 . 



(46) 



(47) 



The parameter space allowed for s e and M3 (or s v ) will 
be determined by doing a simple least square fit. We are 
interested in getting a solution which can lower the %2, 
or 



AX2 =X2(Se,M 3 ) 



xi M 



< 0. 



(48) 



indicating a better fit than the SM. However, for M3 > 
200 GeV, our numerical search did not find any param- 
eter space which can improve the global fitting listed in 
Table [I] This is easy to understand. For Mz s Mz, 
s n <C s e and the corrections in the third column of the 
table is not important. One sees that about half of the 
EW observables get the wrong sign corrections. There- 
fore at most we can only make gain and loss balanced. 

Therefore, we demand that the allowed parameter 
space does not make the fitting worse and this is shown as 
the lower band in Fig. [21 The allowed parameter region 



1.0 

10- 1 
10- 2 
10- 3 
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500 



1000 1500 
M 3 (GeV) 



FIG. 2: The bound on s e and M 3 from EWPT. The upper 
band is the excluding region by too large a deviation from 
SM, \2 > 2xf M - The lower band in the parameter space 
gives comparable to SM results in the global fit. And the 
middle one is the allowed region where (A X2 /xf M ) < 0.01. 



is approximately given by 



M 3 > 166|s e |TeV or |s e | < 0.006 



However, if we relax the global 



A/ 3 



(49) 



lTeV / 
fitting a little, 



(Ax2/xf M ) < 0.01, the constraint can be much looser 
(see the middle band in Fig. 0). Hereafter, the ceiling 
boundaries of the middle and lower bands will be referred 
as 1%EWPT and 0%EWPT respectively. 



IV. THE LHC AND THE ILC 

Here we calculate the Z s Drell-Yan processes at the 
LHC given the constraints on its couplings obtained 
above. One needs to fold in the parton distributions of 
u (d) and u (d) inside the proton. This involves QCD cor- 
rections to the parton model, and are extensively studied 
with relatively good theoretical control [l4|- It is also 
very well tested for the Z production at the Tevatron. 
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Assume narrow width approximation, 



M l 



iM z 1/ 



Mz.Tz. 



, (50) 



the expected number of observed events, say by recon- 
structing from the fi + fj,~ X final states, is 



Nz 3 = L<tt(pp — > -^s^ 

~ — Cz Cexp 
s 



-A 



n + irx 
M z 



(51) 



where 



4tt 2 r(Z fl ) 
3 M Zs 
1 



Br(Z s — > ) 



5r(Z s -> mt) + ——Br(Z s -> dd) 

Cud 



(52) 



L is the luminosity, C u d is the ratio of u(d)-parton dis- 
tribution inside the proton. For the pp hadron collider, 
Cud{pp) ~ 2, A(pp) ~ 32, and C(pp) ~ 600 for a very 
wide range of s [l4j . 

In order to select a signal, it is essential to know the 
branching ratios of Z s . The calculation for Z s decays into 
fermion pair is straightforward. For the SM fermions, 
they are (again setting s 2 ^ = 0.2311): 

T(Z S - 



1 9 2 2 



) = —^-Mz. [0.32739^(1-^) 

2l 



T{Z S -> dd) 



T(Z S -» ee) 



-0.129957s,,c^s e + 0.43022s;] , (53) 
= 4-#-^ s [0.09629 S 2(l- s 2) 

Z47T Cpy 

-0.277396s,c„s e + 0.55451s 2 ] , (54) 



1 92_ 

.2 

-vv 



247rc 2_; M ^[°- 28888s e(l-4) 



-0.09293s t? c I) s e + 0.12571s 2 ] , (55) 

1 <£ 

2 

IV" 



— — Mz s [0.05777^(1 -4 

Z47T C-— 



-0.240364s 7) c^s e 



0.25s 



'/J 



(56) 



After crossing the tt threshold, Z s can decay into a pair 
of t-quarks. We find 



If Mz s is heavier than (Mz + M h o ) , and the kinematics 
are favorable, there is a new decay channel opening up, 



T{Z S -> Zfe§) 



gjc 2 
1927rc^ 



(c, f + sw^s,,) 2 ^ - svvs e c ?) ) z 



M| (59) 



xMz s Vl-4z(l + 8z), z 



Ml 



where we have ignored the mass difference between h\ 
and Z to simplify the expression. The branching ratio of 
shadow Z as a function of its mass is displayed in Fig. 
In generating the figure we have treated Mz s and s € as 



Branching ratio (%) 
17.5 















Z, -. e+ e- 






Z, — rf<7 




Z, -. vii 






Z, ->Zhi 



M z , (GeV) 

300 500 700 1000 1500 

FIG. 3: Branching ratio for the Z a decays as functions of 
Mz s ■ The curves shown here are generated with M h o = 120 
GeV and c a , the Higgs mixing angle as defined in Eq. 117L 
set to one. 

independent parameters and used a very small value of 
s e = 10~ 3 . This is consistent with the bound of Fig. [2 

In the large Mz s limit, say > 1.0 TeV, the s 2 -terms 
will dominate and we obtain a very simple expression for 
the decay width: 



2. 37 ^|fi = 0.1742 
247^ 



M z 



I TeV J V 0.01 



GeV. 
(60) 

We see that for such a heavy Z s its width is indeed very 
narrow, and the various branching ratios are approxi- 
mately given by 



T(Z S — + uu) 



width 
T(Z S 





B u 


= B C 


= B t c 


± 13.81%, 


-4/3[l + 0.41176/? + 0.7979/3 2 ] (57) 


B d 


= B S 


= B b c 


t 4.06%, 




B e 


= B» 


= B T c 


t 12.19%, 


) 2 and the /3 2 term comes from the 


Bu e = 


Bu u = 


--B Vr c 


± 2.44%, 


The decay into a W + W~ pair has 




Bw+w- - 


t 1.219%, 








Bhz - 


s 1.219c 2 % 



w + w~ 



192tt M, 4 



x(l + 20y+12y 2 ), y = 



Mi 



(58) 



(61) 



For SM fermions, the branching ratio is roughly propor- 
tional to (Y^+yj 2 .). It's interesting to note that Z s prefers 
to decay into u-type quarks and charged leptons than 
other SM fermions. This is very different from the SM Z 
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decay. If the ILC is available, and the Z s is found, this 
unique prediction may be tested. 

At the LHC with the center of mass energy y/s = 
14 TeV and using the bench mark luminosity of L = 
100 fb _1 , we calculated the expected number of events of 
Z s into different decay modes by simply folding in the 
branching ratio we obtained earlier and taking into ac- 
count of the phase space factors. 

More importantly, one has to carefully take into ac- 
count the maximally allowed s e obtained from the global 
fit of low energy precision measurements as given previ- 
ously. The expected number of events depend on how 
restricted we are in taking the EWPTs. Fig. 0] shows 
the sensitivity to the A%2/xf M which can lead to two 



V SM\ = o v 1 (1-5 




1.2 1.6 

Mz.(TeV) 



FIG. 4: The maximal expected number of Z s events at the 
LHC for integrated luminosity of 100/6 -1 For a fixed Mz s , 
we have used the largest allowed s e comes from the global fit 
studied in previous section. The left and right panes are for 
0%EWPT and 0%EWPT respectively. 

orders of magnitude difference in the signature. Notice 
the dipping of the signals for smaller Mz s . This is due to 
the much smaller values of s e allowed for these relatively 
light Z s . 

One way of distinguishing between different extra Z 
models will be measuring the branching ratios into dif- 
ferent fermion species. The Z s has the feature that it 
has a relatively large branching ratio into charged lep- 
tons and the t-quarks. For sufficiently heavy Z s the two 
branching ratios are almost equal. Whether this can be 
used as a diagnostic tool at the LHC depends on the 
t-jets and c-jets efficiencies. 

The success of LEP has demonstrated that e + e~ col- 
liders are powerful machines for studying neutral gauge 
bosons. Indeed the search of extra Z bosons has been 
conducted at LEP and the results can be found in ref- 
erence pl|. Looking forward to the ILC, the center of 
mass energy of the collider is lower than that of the LHC 
with y/s = 1 TeV. We can still expect to see extra Z's of 
mass below 1 TeV to be produced. On the other hand, 
facilities designed to have higher luminosity and a bench- 
mark integrated luminosity of 500 fb -1 can be expected. 
Furthermore, the underlying processes involve much less 
QCD uncertainties than in hadronic machines, making it 
a cleaner environment for detecting the extra Z bosons. 
Thus, we can anticipate the branching ratios to be ac- 



curately measured. Moreover, the Z s will be too narrow 
for the total width to be measured. However, spikes will 
be seen at the mass where LHC "discovered" the new 
state. In Fig. [S] we display the result of such a hypo- 
thetical occurrence of a Zg of mass Mz„ — 500 GeV 
and s e — 0.066 corresponding to the maximal allowed 
value from the 1%EWPT fit. The familiar SM Z boson 
resonance peaks sit on the left hand side. A new spike 
appears at Mz,- We magnify the event shape around 
y/s = 0.5 TeV and we see the characteristic dip at the 
left base of the peak corresponding to an extra Z. This 
dip is due to the negative contribution from 7 — Z s and 
Z-Z s interference. Although the resonance factor Kz B z a 
dominates over Kz sJ , Kz 3 z around Mz s , the Kz s z s gets 
an extra suppression in couplings compared to Kz s -y 
and Kz a z which makes the dip visible. 

For 0%EWPT fit, not shown, the spike is not as pro- 
nounced and the width is thinner, thus its studies at the 
ILC will be more challenging. Similarly, the ratio 



<7[pb] 





10 2 


f(pb) 






10 






■ " / / 


A 1 


ft 




\m J 












~460 480 500 520 540" 
Vi(GeV) 1 


r M 


z, = 500 GeV 1 




= 0.066 





0.1 0.5 1.0 

sfi (TeV) 

FIG. 5: The cross section for e+e" -* ff with a 500 GeV Z a 
and s £ = 0.066. (In the main frame, the spike tips have been 
chopped to avoid overlapping among curves.) 



R 



had 



Br(e~ 



qq) 



H+fi ) 



(62) 



where q sums over all quarks except the top, has a pro- 
nounced spike for 1%EWPT fit. This is shown in Fig.© 
The inlay magnifies the region around Z s and illustrates 
the expected interference pattern of two spin-1 particles 
is clearly discernible. The unique event shape is char- 
acteristic of this model which may be used to discrimi- 
nate it from other extra Z models. Finally, the forward- 
backward asymmetry of muon v.s. ^/s is shown in Fig.0 



V. CONCLUSIONS 

We have studied in detail a simple model with an extra 
neutral U(l) boson, dubbed shadow Z. The shadow Z 
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FIG. 6: The ratio of e + e — > 99 over e+e ~* 

with a 500 GeV Z s and s e = 0.066. 
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FIG. 7: The forward-backward asymmetry of muon v.s. y/s 
with a 500 GeV Z a and s E = 0.066. 



mixes with the SM U{l)y gauge boson kinetically which 
is parameterized by s e . The Higgs sector also contains a 
scalar field that interacts with the SM Higgs field. There 
is no direct coupling between the shadow Z and the SM 
fcrmions. This simple model can easily be embedded in 
more elaborate model. Since our motivation is purely 
phcnomcnological we leave this aspect to a future study. 
Instead we embark on a detail analysis of the EWPT 
constraints and other low energy precision measurements 
on Mz B and s e . It is well known that s e can vary greatly 
from model to model. We found that the data constrain 



it to be very small for a wide range of extra Z boson 
masses; see Eq. (|^|) . We conclude that the data favors 
those models in which s e is radiatively generated. 

Not surprisingly the production of Z s at the LHC and 
the ILC depends crucially on s e . In order to ascertain 
whether the signals are observable we define a figure of 
merit measure given by (see Eq. (I48|) ) which we 

found to be positive. We conclude that the shadow Z 
does not give a better global fit to the EWPTs. However, 
we use A\2 to quantify the data tolerance to Z s . We 
found that (Ax 2 /xf M ) ~ 0.01 will lead to an observable 
production of Z s via the Drell-Yan process at the LHC. 

To distinguish the shadow Z from other extra Z mod- 
els (see 0) one has to do as many branching ratio mea- 
surements as possible. The shadow Z has almost equal 
branching ratios into u-type quarks and charged leptons. 
It also has a decay channel into the SM like Z and Higgs 
boson although it is only 1.3%. Similarly for the decay 
into W + W~ pairs. For this we find that the ILC will 
be invaluable for pinning down the nature of the extra Z 
boson. 
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Note Added 

After the completion of this paper our attention was 
drawn to an earlier work that had looked at similar Z' 
models [l6|, and also the Stueckelberg Z' model which 
has almost identical collider signatures |17| and a similar 
electroweak fit |18|. We have checked that our results 
agree where they overlap. A variant of the model has 
also been used in a recent leptogenesis study |l9j . 
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